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ABSTRACT 



We present a detailed structure function analysis of the UV variability of Q0957+561. From optical observations in the 2005-2007 
seasons, we constructed normalized structure functions of the quasar luminosity at restframe wavelengths A ~ 2100 A and A ~ 2600 
A. Old optical records (1995-1996 seasons) also allow the structure function to be obtained ^t A ~ 2100 A, but 10 years ago in 
the observer's frame. These three structure functions are then compared to predictions of both simple and relatively sophisticated 
(incorporating two independent variable components) Poissonian models. We do not find clear evidence of a chromatic mechanism 
of variability. From the recent data, ~ 100-d time- symmetric and ~ 170-d time-asymmetric flares are produced at both restframe 
wavelengths. Taking into account measurements of time delays and the existence of an EUV/radio jet, reverberation is probably the 
main mechanism of variability. Thus, two types of EUV/X-ray fluctuations would be generated within or close to the jet and later 
reprocessed by the disc gas in the two emission rings, al A ~ 2100 A and A ~ 2600 A. The ~ 100-d time-symmetric shots are also 
responsible for most of the ~ 2100 A variability detected in the old experiment (10 years ago). However, there is no evidence of 
asymmetric shots in the old UV variability. If reverberation is the involved mechanism of variability, this could mean an intermittent 
production of high-energy asymmetric fluctuations. The old records are also consistent with the presence of very short-lifetime (~ 10 
d) symmetric flares, which may represent additional evidence of time evolution. Despite these exciting findings, we cannot rule out 
the possibility of unfortunate gaps in the old light curves (or a relatively short monitoring period) and very short-timescale systematic 
noise. We also discuss the quasar structure that emerges from the variability scenario. 

Key words, gravitational tensing - black hole physics - quasars: general - quasars: individual: Q0957+561 



1. Introduction 

Structure functions for optical/UV variability of quasars al- 
low investigation of processes that cause intrinsic variations 
, and t h e properties of intrinsic fl a res (e.g. [K awaguchi et al. 



1998t ICid Fernandes et al. ) 120001: IVanden Berk et al. , ,2004t 



de Vries et al . 2005; Wilh ite et al. I l2008l). flie composition of 



inter vening systems (e.g. lLewis & Irwin 111 9961; ISche chter et al. 
2003j), and the nature of the intergalactic medium (e.g. Hawkins 



20021) . The most recent and comprehensive study of ensemble 



structure functions of non-lensed quasars has been presented by 
[Wilhite et al. (2008). In this work, about 2500 quasars at a me- 
dian redshift of ~ 2 were classified into six different groups ac- 
cording to t heir black hole mass and continuum luminosity. For 
each group. IWiUiite et aTl (12008) derived the square root of the 
noise-less structure function for the ugriz bands. They detected 
the well-known anticorrelation between luminosity and variabil- 
ity, as well as a correlation betw een variability and black hole 
mass (see also lWold et al. 1120071) . Taking these findings into ac- 
count, most of the variations at restframe time lags < 1 year 
are probably intrinsic to the quasars. These might be due to 
differences in accretion rate. The analy sis by Li & Cao (2 008) 
also supports this la st claim. However, IWold et al. I (l2007l) and 
IWilhite et al. I (l2008l) point out the absence of a clear correlation 



between variability and black hole mass for restframe time lags 
< 100 d. Therefore, short timescale variations could not be re- 
lated to changes in accretion rate. 

The Wilhite et al. dataset did not allow for a measure- 
ment of the variability of very bright and massive quasars. 
Moreover, the mechanism of variability might differ from lo- 
cal to distant quasars, and other physical properties (apart from 
luminosity and black hole mass) may play a role (e.g., cir- 
cumnuclear activity, presence of a jet. X-ray activity, etc.). 
Thus, detailed structure functions of well-characterised indi- 
vidual quasars are important tools for understanding different 
quasar p opulations. With resp e ct to t he individual sources, for 
example, ICid Fernandes et alri (l2000l) analysed the members of 
a sample of optica lly selected nearby quasars (z < 0.4; see 
iGiveon et al. II 19991) . They used Poissonian models, and derived 
flare lifetimes and other parameters of interest. A Poissonian 
model describes physical scenarios in which the luminosity is 
due to the superposition of a variable component and a constant 
background (the nonflaring part of the emissivity). The variable 
component is made by the superposition of flares, having a given 
profile and occurring at random times iCid Fernandes et alTj 
2000). Cid Fernandes et al. (2000) reported the poor sensitiv- 
ity of their fits to phenomenological models (flare profiles), with 
the exception of exponentially decaying flares. These asymmet- 
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ric flares led to poor fits. ICoUier & Peterson] (1200 Ih also stud- 
ied optical/UV structure functions of 13 local AGN (Seyfert 1 
galaxies) with very good time coverage and resolution. The flare 
lifetimes (using symmetric triangular flares and certain lag inter- 
vals) were r ~ 5-100 days, with the higher mass AGN having 
larger variability timescales. 

Many previous studies focused on the initial logarithmic 
slope of the structure function. For example, using the square 
root of the structure function, a measured initial slope can be 
compared to predicti ons of possible physical mechanisms of 
quasar variability (e.g. KawagucMetal. 1998; Hawkins 2002j). 
The initial logarithmic slope yS should be less than or equal to 
0.5 for square and exponentially decaying flares (Poissonian 
phenomenological models), wherea s symmetric triangular flares 
produce a steep slope [5 > 0.5 ( Cid Fernandes et aTl 120001: 
iGoicoec hea et al. P 2008). One can also discuss simple models of 
the reverberation scenario, i.e., accretion disc flares induced by 
variable EUV/X-ray irradiation. The cellular- automaton model 
is ab l e to reproduce a slop e p ~ 0.4-0.5 dMineshige et al. I 
Il994t iKawaguchi et al. 1 11998). but cellular-automaton simula- 
tions neglect hydrodynamical effects. These might be basic in- 
gredients in the reverberation s cenario. The on e-dimensional 
hydrodynamical simulations by iManmoto et al. I (11996.) led to 
time-symmetric flares, i.e., flares having symmetric rise and de- 
cay. These time-symmetric flares can account for steep growths 
with slopes above 0.5. The starburst mo del (supernova explo- 
sions) produces an initi al slope jS > 0.7 dAretxaga et al. I [19971 
IKawaguchi et al. 1 119981) . and the microlensing model (extrin- 
sic variability caused by collapsed objects passing close to 
the sight line toward a q uasar) leads to a shallow slope /3 ~ 
0.2-0.3 (lHawkinsll2002h . For non-lensed quasars, the slope 
of the square-root noise-less struct ure function over restframe 
lags > 100 d is ro ughly 0.5 (e.g. I Cid Fernandes et al. I lloool: 
IWilhite et all l2008h . On the other hand, local AGN show a 
variety o f initial slopes over restf rame lags < 100 d: y6 ~ 
0.3-0.8 dCollier & Peterson 1 12001|) . For timescales < 60 d, 
ICollier & Peterson I ( 1200 lb also found that optical and UV slopes 
are consistent with each other This supports a common variabil- 
ity mechanism. 

Over the last 10-15 years, several gravitationally lensed 
quasars have been monitored more or less regularly. Some of 
them offer a unique opportunity to study the origin of the in- 
trinsic signal of quasars, since their intrinsic fluctuations are re- 
peated in two or more lensed c omponents with certain time de- 
lays and magnitude offsets (e.g. lSchneider et anil992l) . Optical 
frames {g and r bands) of the double quasar Q0957-H561 
were taken with the Apache Point Observatory (APO) 3.5 
m telescope during the 1995 and 1996 seasons. These pro- 
duced accurate light curve s that basically include int rinsic varia- 
tions (iKundic et al. Ill997h . IKawaguchi et aTl dl998h studied the 
square root of the structure function of Q0957H-561 using the 
APO g-band light curves (photometric magnitudes), finding a 
shallow logarithmic slope of ~ 0.35 over restframe lags < 200 d. 
Unfortunately, these authors did not subtract t he observational 
noise f rom the photometric measur e ments (e.g., Simonetti et al. I 
119851: ICid Fernandes etalj l2000t ICollier & Peterson I 1200 ll) . 
which has a significa nt effect on the me asured variations at the 
shortest lags. iGil-Merino et al. I d2001h found an initial slope 
of the square-root noise-less structure function of Q0957H-561 
in the R band P > 0.5 (timescale < 100 d). Very recently, 
iFohlmeister et al. I d2008l) separated intrinsic from extrinsic vari- 
ations in Q1004H-41 12, and measured a slope of ~ 0.5 over rest- 
frame lags < 1 year (intrinsic str ucture function in the r band). 
From r-band data of O0909+532. lGoicoechea et al. I d2008l) also 



estimated a steep initial slope /3 ~ I over restframe lags < 100 d. 
No simple model was able to accurately reproduce the shape of 
the structure function of the intrinsic lumi nosity of Q0909-H5 32. 

The first lensed quasar Q0957+561 (IWalsh et al. II 19791) is 
probably the best-studied lens system. This has been inves- 
tigated in several spectral regions, including radio, IR, opti- 
cal, UV and X-ray wavelengths. The radio maps of both com- 
pone nts showed the presence of rad io cores an d ~ O'.'l jets 
(Por cas et al. lll98ll:lGorenstein et al. Il l984. 198H lGarrett et al. I 
,19941: iHaarsma et al. I l2001h . iHutchings . (,2003^ reported evi- 
dence for EUV (~ 1100 A) activity along t he radio jets . Apar t 
from EUV emission associated with the jet, iHutchingsl (l2003l) 
also found EUV emission within a radius of 0'.'3, which is as- 
sociated with a circumnuclear environment. Taking into account 
the redshift of the quasar (z = 1.41), optical observations in the 
g and r bands correspond to middle ultraviolet (MUV) emis- 
sion (~ 2100-2600 A). Hence, fluctuations in the observed op- 
tical magnitude trace variations in the UV luminosity of the 
source. Q0957+561 is a very bright and massive object, with 
both AL^iUSOk) and AL^(3000 A) exceeding 10"^^ e rg s"', and a 
blackhole mass of 2-3 x lO** Mq (iPeng et al.l2006l) . T he double 
quasar Q0957-I-561 is also a X-ray bright source (e.g. IChartasI 
,2000ft . but a possible X-ray jet has not yet been resolved. 

Liverpool Quasar Lens Monitoring (LQLM) is a long-term 
project to follow the variability of gra vitationally lensed quasar s 
with the Liverpool robotic telescope (IGoicoechea et al.]|2008l) . 
The first phase of this project (LQLM I) included a monitor- 
ing programme of Q0957+561 in the gr optical bands. The 
new LQLM 1 light curves of Q0957+561 (2005-2007 sea- 
sons) indicated the ab sence of significant extrinsic variability 
(IShalvapin et al. Il2008l) . In addition, using the time delays (be- 
tween quasar components a nd between optical band s) from the 
APO and LQLM I records, Shalvap in et al. I (|2008) suggested 
that most prominent variations in Q0957H-561 are very probably 
due to reverberation in th e accretion dis c around the supermas- 
sive black hole (see also ICoUier 11200 i V In Section 2, the new 
LQLM I dataset is used to accurately trace the shapes of the 
square-root noise-less structure functions of the quasar UV lu- 
minosity. These shapes are closely linked to the nature of the 
UV fluctuations (see here above), so we investigate the mech- 
anism(s) of variability at two different restframe wavelengths: 
A ~ 2100 A (g band) and A ~ 2600 A (r band). By compar- 
ing the new shape with the old (APO dataset) at A ~ 2100 A, 
we also check for time evolution of the mechanism(s) of vari- 
ability. In Section 3, we present our conclusions and discuss 
the quasar structure (temperature profile and source sizes) that 
emerges from the favoured variability scenario. 



2. Structure function analysis 

The APO and LQLM I light curves of the double quasar 
Q0957-H561 do no t show evidence of extrinsic variability 
dKundic et al. Ill997t IShaivapin et al. II200&) . Thus, we directly 
obtain the noise-less structure function of the intrinsic luminos- 
ity at a given restframe wavelength (instead of one structure 
function for each component, a combined re cord and the cor- 
responding structure function are made, e.g., IKawaguchi et al. I 
,1998.) . The combined records that we use in this paper are de- 
picted in Fig. [T] There is no standard form of the structure 
function, but different approaches to the problem. For example, 
while magnitudes, 5' F(ot ), are often used in optical astronomy 
(e.g. IWiMte et al. I l2008l and references therein), monochro- 
matic fluxes or luminosities, SF(F) or S F(L), are more rele- 
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Fig. 1. APO and LQLM I combined light curves. While the 
light curves of Q0957+561B (open circles) are unchanged, the 
Q0957+561A records (filled circles) are shifted by the 417-d 
time delay and optimal magnitude offsets (producing optimal 
overlaps). Top: APO g-band. Middle: LQLM I g-band. Bottom: 
LQLM I r-band. We note that all magnitudes are expressed in 
the SDSS photometric system. 



vant in radio astro nomy or to accurately c ompare with models 
of variability (e.g. ISimonetti el: al. I [TqSSI: ICid Fernandes et al~l 
I2OOOI: ICollier & Peterson II2OOII) . The structure function SF(L) 
at restframe lag At is estimated through the averaged sum 



SF{L) = (1/2A^)^[(10"'''^"'^ - 10""- 



Am,\2 



(1) 



where m are magnitudes, cr = 0.921 x 10 ° "*"'cr, cr are photomet- 
ric uncertainties, and the sum includes A' pairs verifying -t, ~ 
At. Here, L - 10 " are monochromatic l uminosities in conve- 
nient units and t are restframe times (e.g., ICid Fernandes et aTl 
l2000tlGoicoechea et al. II2008I) . This SF(L) describes typical lu- 
minosity variabilities at different restframe lags. We normalize 
the original structure function to the luminosity variance and 
then take the square-root for convenience, i.e., we analyse the 
normalized structure function / = [5F(L)]'^^/(t(L). 

Restframe lags substantially below the restframe duration of 
the records are considered in the analysis. In the initial selection, 
we take At < P/4(l -i- z), where P is the duration of each whole 
combined record (see Fig. [T]i. Later, only lags before reaching 
the asymptotic behaviours (/ < 1) are taken into account. In 
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Fig. 2. Normalized structure functions of Q0957-I-561 from 
Liverpool telescope (LQLM I) records. These structure func- 
tions of the UV luminosities at ~ 2100 A (g band; filled cir- 
cles) and ~ 2600 A (r band; open circles) are accurately de- 
scribed from very short restframe lags (~ 1-10 days) to lags 
when [S F(L)y^^/o-{L) ~ 1 (asymptotic behaviours). The initial 
growths are consistent with slopes /3 ~ 0.6-0.8 (solid lines). 



Fig. 121 using LQLM I observations, the normalized structure 
function at A ~ 2100 A (g band; filled circles) is compared to 
f at A ~ 2600 A (r band; open circles). Both growths seem to 
be consistent with each other With respect to the initial loga- 
rithmic slopes, in Fig.|2]we also show two fits / = AiAr'f with 
= X^/Ndof values close to 1 {Ndof is the number of degrees of 
freedom). The fits over time intervals 3-14 d (g band) and 6-50 
d (r band) lead to/3 = 0.71 ± 0.08 and /3 = 0.63 + 0.02 (Icr inter- 
vals), respectively. These initial slopes disagree with the predic- 
tion of the cellular-automaton disc-instability model, but roughly 
agree with the time-symmetric flares that appear in the one- 
dimensional hydrodynamical simulations (see Introduction). In 
principle, the starburst model could also account for the mea- 
sured slope at ~ 2100 A (g band). The microlensing slope 
j6 ~ 0.2-0.3 ( iHawkinsll2002l) is strongly inconsistent with the 
LQLM I data of Q0957H-561, which is not at all surprising. We 
are studying intrinsic fluctuations, so microlensing does not play 
any role. 

In Fig. [3j the APO (open triangles) and LQLM I (filled cir- 
cles) structure functions of Q0957-I-561 at ~ 2100 A have sig- 
nificantly different initial growths. This may be a consequence 
of evolution in the variability scenario, since both experiments 
(APO and LQLM I) are separated by ~ 1500 d in the quasar 
restframe. In fact, the oldest (APO) brightness records seem 
to incorporate relatively short fluctuations (which would not be 
present in the LQLM I light curves) that would generate the dif- 
ferences between initial growths. The solid lines in Fig.[3]fit the 
APO behaviours in three time intervals (x^ ~ 1). Their slopes 
areyS = 0.81 + 0.12 at At < 5 d,j6 = 0.66 + 0.06 at At < 11 d, 
and /3 = 0.59 + 0.04 at At < 16 d (Icr intervals). 

In order to discuss the mechanism(s) for (M)UV variability, 
and to quantify the spectral and time evolution of that mecha- 
nism(s), we compare the observed shapes of the structure func- 
tions (Figs. |2|-[3]) with predictions of Poissonian models. In a 
first level of complexity (simplest models), we use three phe- 
nomenological models: square flares (SQF), exponentially de- 
caying flares (EDF), and symmetric triangular flares (STF), as 
well as starburst flares (SBF) produced by supernova explosions 
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Table 1. Solutions for the main Poissonian models. 



Model" 


Observed structure function* 






Ti'- (d) 


T2'- (d) 


SQF(1)+STF(2) 


fLQLM /(2100 A) 


1.21 


266 

u.^uu ,,(,20 


32.6 


175.8 




fLQLM ;(2600 A) 


0.65 


170 +""28 


9 +4.2 


122.7 




/^Po(2100 A) 


1.44 


266 

U.Z.UU _o o24 


11 

ll.u ,2 


107.1 


SQF(1)+SBF(2) 


fLQLM ;(2100 A) 


1.19 


706 +"""2 

-0.664 


151.4 !]5|2 


135 +^^> 




fLQLM ;(2600 A) 


0.99 


630 


88.2 


353 




fApo{2m A) 


1.78 


fi t;c +0.042 
-0.046 


99.4 


32 

■-'^ -6 


EDF(1)+STF(2) 


fLQLM ;(2100 A) 


1.25 


0.902 


180.4 


90 6 +""■■' 

-26.1 




fLQLM ;(2600 A) 


0.73 


0.418 


154.8 


97.2 l 




/iPo(2100 A) 


1.80 


n S'>fi +0.032 


135.6 !i« 


10.8 !j-5 


STF(1)+STF(2) 


fLQLM ;(2100 A) 


1.55 


0.204!°;™ 


24.3 


162.3 l\] 




fLQLM /(2600 A) 


0.80 


0.182 


36.0 +5-^ 


123.9 l^'j 




/iPo(2100 A) 


1.31 


270 


10.8 


108.0 +f 3 


STF(1)+SBF(2) 


fLQLM /(2100 A) 


1.18 


0.052 


11 7 

11./ 2.1 


287 !|32 




fLQLM /(2600 A) 


0.72 


0.608 


112 5 ° 


16V 




/^Po(2100 A) 


1.64 


0.646 


114.3 


33 !J 



" SQF = square tlares, EDF = exponentially decaying flares, STF = 
symmetric triangular flares, and SBF = starburst flares. 

* Normalized structure functions (see main text) from Liverpool tele- 
scope (LQLM I) and Apache Point Observatory (APO) data. 
All measurements are Icr intervals. 
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Fig. 3. Normalized structure functions of Q0957+561 from 
Apache Point Observatory (APO) and Liverpool telescope 
(LQLM I) records. The two structure functions of the luminos- 
ity at ~ 2100 A (APO = open triangles, LQLM I - filled circles) 
correspond to two experiments separated by ~ 10 years, i.e., ~ 4 
years in the quasar restframe. The initial logarithmic slopes (see 
solid lines) of the APO trend vary from fi ~ 0.8 over lags < 5 d 
to /3 ~ 0.6 over lags < 16 d. 



dAretxaga etal.lll997t ICid Fernandes et al. I I2OO0I) . For exam- 
ple, EDF is a simple model to describe asymmetric flares, i.e., 
rapid rises and slow declines, or slow rises and rapid declines 
(the structure function cannot distinguish between both vari- 
ants). S TF is a roug h description of the hydrodynamical simula- 
tions by iManmoto et al. I d 19961) and a relatively good approach 
to the X-ray shots found bv lNegoro et al. I (Il994i) . Each model is 
characterised by a shape function s(t), f - [^(t)]' ^^, and these 



shape functions appear in Appendix B o f Cid Fernandes et alTI 
(2000) and Eq. (17) of lAretxaga et al. I (|1997|). The 4 level-1 
models of / only include one free parameter: flare lifetime t 
(e.g., T - T for a square flare of duration T and r - 5tsg for 



a supernova explosion, where fjg is the time when the super- 
nova remnant reaches the maximum of its radiative phase). In the 
Poissonian framework, we can consider more complex schemes. 
Thus, in a second level of complexity, the luminosity is assumed 
to be due to the superposition of a constant background and two 
independent variable components. It can easily be shown that 
/= [wiii(Ti) -t- (1 - wi).S2(t2)]'^^, where wi is the ratio between 
the variance of the first fluctuating component and the total vari- 
ance, and ii and S2 are the shape functions of the two indepen- 
dent components. 

The possibility that AGN variability is caused by several 
independent processes was suggested in previous work (e.g. 
iKawaguchi et aLl 119981: IComer & Peterson I l2001h . Here, us- 
ing 10 level-2 models of / (SQF-i-SQF SQF-i-EDF SQF-nSTF 
SQF+SBF, EDF+EDF EDF+STF, EDF+SBF STF+STF 
STF+SBF, and SBF-i-SBF), we also explore this possibility. 
Each of these 10 possible combinations is characterised by three 
free parameters: wi, and two lifetimes ti and T2 (see above). To 
evaluate the quality of the fits obtained with the full set of 14 
models, we analyse the (reduced chi-square) values. For an 
acceptable fit,;^f^ is expected to be in the range Ndofi^ilNdo/y^^ 
(allowing - Njof differences of up to two standard deviations 
of thex^ distribution), which implies 0.6 < < 1.4 (Njof ~ 45). 
As expected from the measured slopes, the five best solutions in 
terms of x^ are associated with (level-2) models incorporating 
either STF, or SBF, or both of them. However, only one out of 
our 42 fits to /lqlm /(2100 A), /lqlm /(2600 A), and /ato(2100 
A) (3 observed shapes x 14 models) gives < 1.4. This result 
means that the uncertainties are slightly underestimated or alter- 
natively, the models do not describe all details. The statistical 
uncertainties in each structure function, see Eq. (1), are com- 
puted as the standard deviations of the means (averaged sums) 
for the different time lag bins. Although we use a ve ry popular 
estimator of uncertainties. ICollier & Peterson I (1200 ll) noted that 
not all pairs of data in a given bin are independent. To address 
this problem, they multiplied their errors by 2''^. We adopt the 
Collier & Peterson's perspective and fit again the observed struc- 
ture functions. 



Goicoechea et al.: Structure function of the UV variability of Q0957+561 



5 



After slightly enlarging the error bars, we concentrate on 
level-2 models consisting of STF/SBF and anything else. Among 
these 7 models of /, we select those giving < 2 for the 
three observed structure functions (main Poissonian models). 
The set of solutions is presented in Table [1] With regard to 
the spectral behaviour, four out of the 5 models in Table [T] are 
able to accurately and simultaneously describe /lqlm /(2100 A) 
and fiQLM /(2600 A), i.e., they fit both LQLM I shapes ac- 
ceptably well (0.6 < < 1.4). The solutions for the first 
model (SQF+STF) indicate that symmetric triangular flares at 
~ 2100 A have a Ufetime longer than the duration of STF 
at ~ 2600 A. However, this is puzzling with the accretion 
disc paradigm. Similar lifetimes are expected in a reverbera- 
tion scenario. Moreover, t(2100 A) < t(2600 A) in a disc local- 
instability scenario, e.g., the thermal timescale is proportional to 
{R^IGMyi^, where M is the mass of the central black hole and 
R oc A^l^ is the emission radius (e.g.. iGoicoechea et al. Il2008l 
and references therein). The solutions for the second model 
(SQFh-SBF) incorporate supernova explosions with fjj(2100 A) 
< fsj(2600 A). However, the timescale tsg is exlusively related 
to the circumst ellar density and the to tal energy released in each 
explosion (e.g.. lAretxaga et alTII 19971) . so a chromatic timescale 
is not the expected result. 

The solutions for the EDF+STF and STF-i-SBF models are 
the most interesting ones (LQLM I data). There is a clear de- 
generacy between EDF and SBF, which seem to work in a sim- 
ilar way. Hence, the solutions for the third and fifth models 
are interpreted as evidence in favour of the coexistence of ~ 
100-d time-symmetric flares and longer time-asymmetric shots 
(see Table [1]). Both kinds of flares at both wavelengths can be 
mostly due to reverberation, i.e., two types of EUV/X-ray varia- 
tion in the vicinity of the accretion disc axis that are reprocessed 
by two annuli of the disc g as. Measurements of time d elays 
dKundic et al. Il 1 9971: ICollier 1 ^0 1 : Shalvapin et al. II2008I) . and 



the existen ce of an EUV j et (iHutchings ,,2003i) and a bright X- 
ray source jChartas Il2000h . support this physical scenario. From 



Table [T] 



168 d and t. 



105 d is a good compro- 



mise between the results at the two wavelengths and for the two 
flare asymmetric profiles (exponentially decaying and starburst 
profile). Although the relative variances depend on wavelength, 
there is no need to invoke some mechanism other than rever- 
beration. While the ~ 2100 A variability is consistent with the 
asymmetric flares producing most of the variance at this wave- 
length (Wflivm - 12.%'), the variance at ~ 2600 A would be mainly 
due to the symmetric flares (wasym - 40%). This chromaticity in 
the relative variances is associated with the difference between 
time coverages, gaps and artifacts in the g (~ 2100 A) and r (~ 
2600 A) bands (see middle and bottom panels in Fig. [TJ. The 
LQLM I g-band combined record lasts longer than the LQLM I 
combined curve in the r band, so it contains a prominent decline 
that is not present in the shorter record. Besides the longer time 
coverage, the g-band curve fills a gap corresponding to the peak 
of an important event. The presence of a few artifacts could also 
play a role. In Fig.|4] the LQLM I shapes are compared with the 
adopted solutions: 'Wasym — 0.72, Jasym - 168 d, and t,,,,,, - 105 d 
for fLQLM /(2100 A), and Wasxm = 0.40, Tasxm = 168 d, and t,„„ 
= 105dfor/LeiM/(2600A). 

In spite of the enlargement of error bars and the use of 
relatively sophisticated modelling, no model leads to an ac- 
ceptable global solution at both wavelengths and both epochs. 
The APO data are only consistent, in terms of x^, with the 
STF+STF model. The corresponding solution includes ~ 100- 
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Fig. 4. Observed structure functions and adopted solutions. For 
the LQLM I experiment, we show EDF+STF laws with param- 
eters: wi = 0.72, Ti = 168 d, and T2 = 105 d (g band; solid line), 
and wi = 0.40, t\ - 168 d, and T2 - 105 d (r band; dashed line). 
SBF-hSTF laws with the above parameters are also able to re- 
produce the LQLM I shapes. For the APO data in the g band, we 
display a STF+STF law with parameters: w\ - 0.27, ti = 1 1 d, 
and T2 = 105 d (dashed-dotted line). 



d time-symmetric flares and other family of very short-lifetime 
(t ~ 10 d) symmetric flares. This last kind of fluctuations is re- 
sponsible for about 1/4 of the total variance. The dominant ~ 
100-d shots are also detected in the LQLM I experiment, so the 
main source of these flares do not evolve over decades (in the ob- 
server's frame). Nevertheless, the absence of longer asymmetric 
flares (APO data) is a mystery. This could be due to unfortunate 
gaps in the light curves, time coverage of the monitoring, etc, 
or simply reflects the evolution of the central engine. Assuming 
that UV asymmetric flares are triggered by EUV/X-ray fluctua- 
tions (reverberation in the accretion disc; see above), we would 
be dealing with intermittency in the generation of high-energy 
asymmetric fluctuations (within or close to the jet). As we com- 
mented here above, the structure function /afo(2100 A) agrees 
with the presence of very short-timescale shots, which are not 
detected from the new LQLM I records. These fluctuations may 
be caused by some kind of observational (systematic) noise or 
alternatively, they might be related to an episode of very short- 
timescale activity inside the disc, the jet or other region. In Fig.|4] 
the adopted solution (STF+STF model with parameters wi = 
0.27, Ti = 11 d, and T2 = 105 d) fits /apo(2100 A) reasonably 
well. 



3. Conclusions and discussion 

We present a novel and rigorous analysis of the structure func- 
tion of the UV variability of the gravitationally lensed quasar 
Q095 7+561. New Liverpoo l telescope data (2005-2007 sea- 
sons; [Shaly^inetair||200l) allow us to construct normalized 
structure functions of the quasar luminosity at two restframe 
wavelengths: X ~ 2100 A (g band data) and X ~ 2600 A (r 
band data). Old Apac he Point Observator y records in the g band 
(1995-1996 seasons; [Kundic et al. 119971) are also used to check 
the possible time evolution of the variability at ~ 2100 A. The 
observed shapes of the structure functions are compared to pre- 
dictions of a large set of Poissonian models. This set of models 
includes the simplest and well-known ones, consisting of only 
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one variable component dCid Fernandes et aT7ll2OO0l and refer- 
ences therein), as well as hybrid models incorporating two inde- 
pendent variable components. 

Several hybrid (or level-2) models are able to account for 
both Liverpool telescope structure functions (see Table[T]). Some 
of them contain flares with unrealistic profile (square flares) and 
lead to solutions that are difficult to interpret. Fortunately, we 
also find reasonable solutions in which ~ 100-d time-symmetric 
and ~ 170-d time-asymmetric flares are produced at both rest- 
frame wavelengths. Exponentially decaying and starburst flares 
(and very probably other time-asymmetric shots) work in a sim- 
ilar way. Therefore, the good behaviour of the starburst ingredi- 
ent does not necessarily implies the existence of supernova ex- 
plosions, but the production of highly asymmetric shots. What 
about the mechanism of intrinsic variability?. The old and very 
recent gr light curves of Q0957 -1-561 led to time delays be- 
tween quasar components and be tween optical bands that mainly 
suppo rt a reverberation scenario (ICoUier 11200 UlShalvapin et al. I 
l2008h . Thus, reverberation would be the main mechanism of 



variability. The presence of an EUV/radio jet (e.g., Garrett et al^ 
1994; Hutchings 2003) and a bright X-ray source (Chartas 
2000 ) also suggests the viability of this mechanism: two types 
of EUV/X-ray fluctuations that are generated within or close to 
the jet, and later reprocessed by two rings of the disc (each ring 
corresponds to a different restframe wavelength). On the other 
hand, one can also justify both kinds of flare profile. For ex- 
ampl e, the cellular-automato n model produces asymmetric shots 
(e.g.. IKawaguchi et al. Ill998i). and hvdrod vnamical simulations 
lead to symmetric flares dManmoto et al. I ll996). 

The ~ 100-d time-symmetric shots seem to be also responsi- 
ble for most of the ~ 2100 A variability detected in the Apache 
Point Observatory experiment, but there is no evidence of asym- 
metric shots in the old UV variability. This absence of asymmet- 
ric flares may be due to gaps in the light curves, a relatively short 
monitoring period, etc. Alternatively, it could means an evolu- 
tion of the central engine, i.e., intermittent production of high- 
energy asymmetric fluctuations. The Apache Point Observatory 
structure function is also consistent with the presence of very 
short-lifetime (~ 10 d) symmetric flares. This kind of flare might 
be caused by observational systematic noise, or perhaps, repre- 
sent additional evidence for time evolution. Our results do not 
support a previous claim for the possi ble starburst origin of some 
events in the old g-band light curves dUllan et al. Il2003h . Despite 
the presence of two twin events (one in each quas ar compo- 
nent) with an anomalous delay dGoicoechea I l2002h . the asso- 
ciated shot probably occurred at the base of the jet or in the 
circumnuclear region, but it was not originated by a supernova 
explosion. 

Very recently, several studies have showed evidence that 
optical/UV variability of quasars on restframe timescales 
> 100 d is mainly d r iven by variations in accretion rate 



20071: lArevaloeta l. 2008: Li & Caol l2008t 



IWiMteet al.]|2008h . Here we are discussing the UV variability 
of Q0957-I-561 at restframe lags < 100 d. Q0957-H561 is a very 
bright and massive object, a nd this population could not be stud- 
ied bv' WilhiteetaTltBOOSl) . The first lensed quasar has also an 
EUV/radio jet (and important X-ray activity; see above), so high- 
energy variations in the surroundings of the disc axis and their 
reverberation are possible. For example, on timescales below 
100 days, the optical variations of the local Seyfert galaxy NGC 
5548 are related to its X-ray variations dCzernv et al. 1 119991) . 
Hence, part of the optical variability of this AGN (timescales < 
100 days) could be explained by X-ray reprocessing. Chromatic 
delays for the local Seyfert galaxy NGC 7469 also reveal a re- 



verberation scenario dColher et al. lfT999l) . lArivalo et al. I (l2008h 
reported on an illustrative example of mixed variability. The lo- 
cal quasar MR 2251-178 has been monitored simultaneously 
in X-rays and optical bands. All spectral regions were signif- 
icantly variable, and the fluctuations were clearly correlated. 
lArevalo et afl d2008l) indicated that pure reprocessing of X-rays 
cannot account for both ~ 100-d and ~ 500-d timescale optical 
variability. They claimed that two distinct mechanisms produce 
the variability: accretion rate variations plus reverberation, and 
the shortest timescale optical events are due to reverberation. 

For an irradiated disc, in general, we obtain a temper- 
ature profile shallower th an the standard one T oc r^^/'^ 
(IShakura& Sunvaev ll l973l) . The reverberation hypothesis as- 
sumes that the optical/UV disc regions are irradiated by EUV/X- 
ray photons from the vicinity of the disc axis. If the high-energy 
source is placed on the axis and at a height Hx above the thin disc 
(disc thickness << Hx), the non-standard temperature profile is 



T(r) 



3GMM {\-A)LxHx 



87rcrr3 47rcr(//2 + r^fl^ 



1/4 



(2) 



where G is the gravitation constant, cr is the Stefan constant, 
M is the mass accretion rate, A is the disc albedo, i.e., the ra- 
tio of reflection to incident high-energy ra diation, and is the 
luminosity of the irradiating source (e.g.. ICackett et aini2007l 
and references therein). Moreover, considering r » Hx (and 
thus, a standard temperature profile), the typical radius of the in- 
tensity distribution at a given restframe wavelength A should be 
greater than the standa rd value (for standard structure, see, e.g., 
IShalvapin et"ani2002l) . Eq. (2) leads to 
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(1 -A)LxHx 


1/3 


kA 
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where k is the Boltzmann constant, h is the Planck constant, 
and c is the speed of light. This non-standard typical radius is 
produced by both the heating due to irradiation and the viscous 
heating in the disc. Finally, we point out that shallow tempera- 
ture profiles (from reverberation) could be co nsistent with mi- 
crolensing data of some lensed quasars (e.g., iPoindexter et al. I 
2008). Moreover, the n on-standard sizes of several lensed and 
microlensed quasars (Mor gan et al. Il2007t iPoolev et al. Il2007h 
might be related to relatively high irradiation-to- viscosity ratios 
IV R = 2(1 -A)LxHx/3GMM. 
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